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Stability issues in the modified Starobinsky model 
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We discuss the stability of the anomaly-induced inflation (modified Starobinsky model) with respect to the 
arbitrary choice of initial data and with respect to the small perturbations of the conformal factor and tensor 
modes of the metric in the later period of inflation and, partially, in the present Universe. 
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The basic principles of the anomaly- induced 
inflation has been explained in [[I] (see also refe- 
rences and notations therein). 

The advantage of this inflationary model is that 
it requires smaller amount of the cosmological 
phenomenology than the usual inflaton models. 
To some extent, this model is mainly based on 
the principles of quantum field theory. In par- 
ticlar, the simplest inflationary solution follows 
from the anomaly-induced quantum correction to 



the vacuum action 
the equation 
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which is a solution of 



f = (wC 2 + bE + cDR) . (1) 
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The coefficients w, b, c are defined in [0, C 2 is 
a square of the Weyl tensor and E is an integrand 
of the Gauss-Bonnet topological term. This equa- 
tion admits an explicit solution [ |2| . It is easy to 
see that this solution contains an ambiguity be- 
cause an arbitrary conformal functional Sclg^u] 
plays the role of the "integration constant" for 
the Eq. (0. However, this "integration constant" 
does not affect the equation for the conformal fac- 
tor of the metric and in this respect the initial in- 
flationary solution follows from the exact effective 
action. The tempered form of expansion which is 
observed at the later inflationary phase is not ex- 
act, but it has quite a robust background. 

The stability of the inflationary solution from 
the initial stage until the graceful exit and the 
stability of the classical solution in the theory 
with loop corrections represent a strong consis- 
tency test of the model. Let us start from the 



initial stage of inflation, when the particle con- 
tent Nq, Ni/2 N\ (number of scalar, fermion and 
vector fields) of the theory provides the stability 
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of the exponential solution 
a(t) = a ■ exp^Hgt) 
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of Starobinsky [Ej. According to [0], the condi- 
tion ((2J) is independent on the cosmological con- 
stant and on the choice of the metric k = or 
k = ±1. For the sake of simplicity we consider 
k = and also assume that the cosmological con- 
stant is small during inflation, such that the last 
is driven by the quantum effects only. 

The original Starobinsky model deals with the 
unstable case. The initial data are chosen very 
close to the exponential solution © such that the 
inflation lasts long enough. Using the 0-0 com- 
ponent of the Einstein equations with quantum 
correction, Starobinsky constructed the phase di- 
agram of the theory. This phase diagram repre- 
sent several distinct attractors, FRW behaviour is 
one of them and others represent physically un- 
acceptable run-away type solutions. In the mod- 
ified version of the model [0, the inflation starts 
in the stable phase J2J). In this case the phase 
diagram, dual to the one of [ has the form 
shown at the Figurel. This phase portrait indi- 
cates to the unique stable solution . Therefore 
the anomaly-induced inflation does not depend 
on the choice of initial data and the possible run- 
away type solutions do not represent a threat. 
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Figure 1. The phase diagram for the stable ver- 
sion of the Starobinsky model, corresponding to 
the MSSM particle content. All the variables are 
the same as in [Pf. 



Let us now consider another extreme case, that 
is the present Universe dominated by the positive 
cosmological constant A. Then the equation of 
motion for the conformal factor a(t) has another 
approximate solution 

o(t) = a Q • exp(iJ s i) , H c = (4) 

The present cosmic scale is defined by the magni- 
tude of the Hubble parameter H sa 10 -42 GeV. 
Then the unique active degree of freedom which 
contributes to the vacuum quantum effects is a 
photon, and hence (No, Ni/ 2 N\) = (0, 0, 1), such 
that the inflation is unstable. But this is not all 
the story, because the stability of the low-energy 
solution does not follow from the instability 
of the high-energy solution @. Let us remem- 
ber that the corresponding equations of motion 
(see, e.g. [ 0J Q) include the fourth derivative 
terms and therefore the stability of the solution 
Q with respect to the small perturbations of the 
conformal factor a(t) — lna(i) can not be seen 
as a trivial occurrence. It may happen that some 
run-away type solution is stable instead. Hence, 
at this point we meet a very strong consistency 
test for the whole approach. 

Consider the perturbation H — ► H + A cxp(Xt) 
where A = const. The solutions of the character- 



istic equation for A have the form [0] 

Ai = -4tf, \ m = - i H ± ^=i. (5) 
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It is nice to see that the present-day Universe is 
not in danger due to the higher-derivative terms, 
because the positive cosmological constant pro- 
vides stability of the solution Q . 

The next test is related to the stability of the 
inflationary solution at the last phase, when the 
quantum effects of massive fields temper the ex- 
ponential behavior. In this case we can use the 
approximate analytic method or numerical simu- 
lations. The results of both methods are the same 
[0]. Let us briefly describe the analytic method. 
The stability or instability with respect to the 
small perturbations depends on the behavior of 
a(t) at the relatively small intervals of time, when 
the Hubble parameter H can be treated as a con- 
stant. Of course, when we move from one such in- 
terval to another, H changes providing a source 
for the perturbations. The direct calculations 
give the following equation for the perturbations 
a — > cr + y(r), where we used "renormalized" time 
variable r = t/H, H = const: 

y+7'y + 2(6- ~)y- — y- — fy = (6) 

\ cJ c c 

At this point we assumed, as before, a relatively 
small value of the cosmological constant. The 
last equation has a very special form, because all 
the coefficients are constants and all but the last 
have the magnitude of the order one. The last 
coefficient is extremely small because of the fac- 
tor / <C 10~ 9 . The stability of equation with 
constant coefficients may be explored, e.g. using 
the Routh-Hurwitz (RH) conditions. A priory the 
RH determinants may have an arbitrary sign, but 
in our case they all turn out to be positive, such 
that the stability of the tempered inflation holds 
until the region H as M*. 

Let us consider the stability with respect to 
tensor perturbations of the metric. In the co- 
variant formalism (see, e.g., the evolution of 
the tensor degree of freedom is described by the 
coordinate-dependent scalar factor h(t, r) of the 
tensor mode. The dynamical equation for h(t) 
is very complicated [ El EI j even for the theory 
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of a massless fields. Moreover, this equation con- 
tains an ambiguity due to the conformal func- 
tional Sclg^u], which we discussed above. One 
can fix this ambiguity by choosing the proper vac- 
uum for the perturbations [ [7] • As a result we 
meet an almost fiat spectrum of the perturba- 
tions, however their amplitude may increase very 
fast. At the same time the amplification of the 
amplitudes performs slower than the expansion 
of the conformal factor. As a result the total me- 
tric becomes more and more homogeneous and 
isotropic. This is a situation at the initial stage 
of inflation. 

At the last stage, e.g. in the last 65 e-folds, the 
equation for h(t) is greatly simplified due to the 
enormous number of the total e-folds between the 
beginning and the end of inflation. The typical 
value for a depends on the model (see [H]) but it 
varies between 10 10 and fO 32 . Obviously, the a 
itself may be treated as a very large number in the 
last 65 inflationary e-folds. This feature greatly 
simplifies the equation for h(t,r). In particular: 

i) The conformal functional S c [g^ v ] and the 
choice of the classical action of vacuum have 
no importance. The equation for hit, r) is 
completely defined by the universal /3-functions 
w, b, c for the vacuum parameters. In particular, 
the difference between the equations of [Ej and 
[Ej (it is due to the different choice of <S C [<?,«,]) 
disappears in this approximation. 

ii) The terms with space derivatives Wh(t, r) 
are suppressed by the factors of cxp(— 2a) and 
therefore are negligible 

iii) The remaining equation for hit) has the 
form (in the same time variables as in ©) 

"h +6 h +11 h +6 h - — h = , (7) 

W<7f 

where af is a value of conformal factor af « 
a(tf) = 1/f corresponding to the point of tran- 
sition from stable to unstable inflation. Since the 
magnitude of a is huge, one can disregard its 
variation in the last 65 e-folds. 

It is remarkable that the general structure of 
the CEQ. (0 is quite similar to the one of the 
Eq. © for the perturbations of a{t). But it 
is even more remarkable that the solution of 



does not have growing modes. One of the roots 
of the characteristic equation is of the order of / 
and others of the order of one, but all of them 
have negative real parts. For the physically rea- 
sonable choice of the initial data the amplitude of 
the perturbations almost remains constant. 

Finally, we can conclude that the stability of 
the inflationary solution with respect to the per- 
turbations of the conformal factor and the tensor 
mode of the metric holds from the initial stage 
(when the quantum fields may be approximately 
considered massless) until the scale M* , when the 
most of the sparticles decouple and the inflation 
becomes unstable. Indeed, the stability of the 
model may be jeopardized in the transition pe- 
riod, where we expect to meet rapid oscillations of 
the conformal factor which should lead to reheat- 
ing. The main line in the further development of 
the model must be related with the quantitative 
description of the decoupling and the transition 
period, investigation of a reheating and density 
perturbations. Indeed, the last issue has been 
already considered [ |H] in the framework of the 
original Starobinsky model [[3]. 
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